Large mass dileptons from the passage of jets through quark gluon plasma 
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We calculate the emission of large mass dileptons originating from the annihilation of quark jets 
passing through quark gluon plasma. Considering central collisions of heavy nuclei at SPS, RHIC 
and LHC energies, we find that the yield due to the jet-plasma interaction gets progressively larger 
as the collision energy increases. We find it to be negligible at SPS energies, of the order of the 
Drell-Yan contribution and much larger than the normal thermal yield at RHIC energies and up to 
a factor of ten larger than the Drell-Yan contribution at LHC energies. An observation of this new 
dilepton source would confirm the occurrence of jet-plasma interactions and of conditions suitable 
for jet-quenching to take place. 



I. INTRODUCTION 

Collisions of relativistic heavy ions are performed with 
a hope to create a plasma of quarks and gluons (QGP) 
and to study the deconfined state of strongly interacting 
matter. Theoretical investigations of the formation, evo- 
lution and signatures of the QGP have played an impor- 
tant role in planning the experiments for these studies. 
One such investigation concerns the radiation of dilep- 
tons which are considered as penetrating probes of such 
matter due to their very long mean free path. Several 
sources of large mass dileptons have been identified and 
studied. These include the radiation of Drell-Yan pairs 
from the annihilation of primary partons in the nuclei [Q , 
thermal dileptons from the annihilation of quarks and an- 
tiquarks in the plasma and from hadronic reactions 
in the hot hadronic gas^| following the hadronization 
of the plasma. Dileptons from the correlated decay of 
charm and bottom quarks have also been considered [Q. 
Dileptons from the decay of hadrons like p and will pro- 
vide information about the possible medium modification 
of hadrons. The pre-equilibrium production of dileptons 
has also been investigated by several authors in various 
approximations |||, ||, 0, ||, |9| . 

Experimentally the low mass dilepton excess observed 
in the nucleus-nucleus collisions at SPS energies [0 has 
lead to suggestions that interactions among hadrons, and 
perhaps their properties as well, are altered in hot and 
dense nuclear matter The interpretation of the ex- 
cess observed in intermediate mass dileptons has 
been done in terms of a thermal hadronic environment 
with some quark-gluon plasma component [ij, p^ . 
The measurement of dileptons has already made a con- 
siderable progress at RHIC experiments . 

In the present work, we discuss the emission of large 
mass dileptons induced by the passage of high energy 
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quark-jets through the QGP. The general formulation 
for their emission rate along with with a discussion of 
the evolution of the plasma, including the possible chem- 
ical equilibration, and the spectrum of jets and Drell-Yan 
dileptons, is given in Sect. II. The numerical results and 
a discussion are given in Sect. Ill followed by a summary. 



II. FORMULATION 

Let us consider a collision of two heavy nuclei at rela- 
tivistic energies. A large number of minijets produced in 
the collision and the subsequent rescattering of the par- 
tons will lead to the formation of a thermalized plasma 
at some initial temperature Tq at time tq. This plasma 
may or may not be in chemical equilibrium. 

The high energy partons, which are produced at times 
T ~ 1/pt Tq, will traverse through this plasma, losing 
energy through collisions and radiation of gluons, a phe- 
nomenon widely known as jet-quenchin g a nd considered 
to be an important probe of the QGP [0 . 

The quark jets passing through the plasma will pro- 
duce large mass dileptons by annihilation with the ther- 
mal antiquarks. They will also produce high energy pho- 
tons ||l8| . If these additional dileptons and photons are 
detected, they will provide a confirmation for the forma- 
tion of the plasma and the occurrence of jet-plasma inter- 
actions and prepare the ground for a deeper understand- 
ing of the jet-quenching phenomenon and the evolution 
of the colliding system. We shall see that the produc- 
tion is quite sensitive to the temperatures at early times 
and controlled by the initial distribution of high energy 
partons. 

In order to proceed we assume that the phase space 
distribution of the quarks in the medium can be decom- 
posed as 

/(p)=/th(p)+/jct(p) (1) 

where we have a thermal component that is characterized 
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by a temperature T, 



/th(p) = cM-E/T). 



(2) 



and a jet component which is assumed to be limited to 
partons having a ^ 1 GeV. This separation is kine- 
matically reasonable, as the jet spectra fall off as a power 
law and can thus easily be differentiated from thermal 
components. We shall, somewhat arbitrarily, limit them 
to having pt > 4 GeV. In addition we give some results 
when only jets with pt > 6 GeV are considered. 

The phase-space distribution for the quark-jets is given 
by the perturbative QCD result for the jet yield as 



/jct(p) = 



1 



gq ttR^ tpt coshy (Pprdy 
X S{r, - y) e(T - n) e(w - t) Q{Ri_ - r) 



(3) 



where (/^ = 2 x 3 is the spin and colour degeneracy of the 
quarks, i?x is the transverse dimension of the system, 
Ti ~ 1/pt is the formation time for the jet and r] is the 
space time rapidity. We shall take T„iax as the smaller 
of the lifetime of the QGP [vf) and the time taken by 
the jet produced at position r with a velocity v in the 
transverse direction to reach the surface of the plasma. 
The distance covered by the parton during this passage 
is given by 



(4) 



where cos = v • r . We assume here that the parton is 
massless and travels with the speed of light. The boost 



J 



invariant Bjorken correlation pO| | between the rapidity y 
and 77 is assumed. 

These phase space distributions can be folded with the 
cross section of the annihilation of the quarks to obtain 
the rate of emission of the dilcptons. The rates are then 
convoluted with the space time history of the evolution 
of the plasma, using a standard procedure. 



A. Rate of emission of dileptons 

The cross section for quark-antiquark annihilation into 
dileptons {ah ^ t^t~ ) is given by 
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T M2 



(5) 



where the sum runs over the flavour of quarks, Nc — 3, 
and s and e/ stand for the spin and the charge of the 
quark. 

Using kinetic theory the reaction rate for the above 
process can be written as 



R 



d^Pa 



fa{Pa) 



fb{Pb)aiM^)v,,, , (6) 



(27r)3-'"'^"'y 27r3 
where fi stands for the phase-space distribution of the 
quark or the antiquark, and pb are their momenta 
respectively and the relative velocity is 



EaEb — PaPb 
EaEb 



(7) 



After some algebra pTl this can be rewritten as 



dR _ a{M^ 
dM^ ~ ~ (27r)6 



Xa dXa d(j)a Xb dxb d(l)b dya dybfa fb 6 [M^ - 2M^XaXb COsh(ya - yb) + 2M'^XaXb COS (j)b] (8) 



r 



where Xa = p^/M, Xb — p^/M and ya and yb are the 
rapidities. The integrations over the azimuthal angles 
yield 



dR 

Zm2 



such that, 



(27r)5 



Xa dXa Xb dxb dya dyb fa fb 

-1/2 



4.XaXb - {2a;aa;(,cosh(?/a - yb) - 1}' 



(9) 



1 < 



2XaXb C0Sh(2/a - 



- 1 



< Xa^b < 00 , 

-00 < ya,b < 00 . 



< 1 



(10) 



When fa and fb are given by a thermal distribution 

/th(p) - exp(-S/r) = exp(-pT cosh y/T) (11) 
the above integral can be performed |Q to obtain 

dR _ a{A'P) 



dAP 2(27r)^ 



-APTKiiM/T) 



(12) 



The expression in Eq. (^) can be used for phase space 
distributions of arbitrary form. When the distributions 
fa and fb are taken as thermal, we shall call the yield 
as thermal and when one of them is thermal and the 
other one is from a jet, we call the yield as due to the 
passage of jets through the plasma in an obvious nomen- 
clature. Noting that the jet-distribution (Eq. ^) depends 
on the quark flavour, we take charge- weighted average of 
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To 


To 




A, 






(fnr/c) 


(GeV) 






(GeV/fm^) 


Bjorken formula 


SPS 


0.20 


0.345 


1.0 


1.0 


26.10 




0.50 


0.254 


1.0 


1.0 


7.94 


RHIC 


0.15 


0.447 


1.0 


1.0 


72.38 




0.50 


0.297 


1.0 


1.0 


14.43 


LHC 


0.073 


0.897 


1.0 


1.0 


1172 




0.50 


0.473 


1.0 


1.0 


90.92 


Self-screened parton cascade 


RHIC 


0.25 


0.67 


0.34 


0.064 


61.4 


LHC 


0.25 


1.02 


0.43 


0.082 


425 



TABLE I: Initial conditions for the hydrodynamical expan- 
sion estimated from multiplicity densities using the Bjorken 
formula and from a self-screened parton cascade model |2q]. 



the distributions for u, d, and s quark-jets and for the 
corresponding anti-quark jets, while using /jet in Eq. (|^): 



We take the particle rapidity density as 750 for central 
collisions of lead nuclei at SPS energies. At RHIC ener- 
gies we estimate the particle rapidity density as ~ 1260, 
based on the measured charged particle pseudo-rapidity 
density for central collisions of gold nuclei at ^/snn — 
200 GeV. For central collisions of lead nuclei at LHC en- 
ergies we use dN/dy = 5625, as suggested in pl| |. 

The initial temperatures and times from these esti- 
mates are given in Table | under the heading Bjorken 
formula. We also use an alternative estimate where the 
formation time is taken as 0.50 fm/c, so that the corre- 
sponding initial temperatures are considerably reduced. 
We shall see that this drastically alters the yield of dilep- 
tons from interactions within the plasma, while the yield 
due to the jet-plasma interaction only reduces by a factor 
of about 2. 

We further approximate the nuclei as being spherical in 
shape with a uniform density and assume that the energy 
density at the transverse position r is proportional to the 
number of collisions so that 



e(To,r)cx4p2 [rI ~ 



(15) 



f 



;p)/Ve^ 



(]^3^ where po is the density of the nucleons. The normaliza- 
tion is then determined by the condition that 



where the sum runs over the flavours of quarks (anti- 
quarks). The thermal distributions are assumed to be 
flavour- independent, so that the relation (Eq. |ll|) is re- 
covered for them. 



B. Evolution of the plasma 

The radiation of dileptons from an expanding and cool- 
ing plasma including the contributions from a mixed 
phase and the hadronic phase with transverse expansion 
of the plasma has been studied in great detail (see e.g. 
Ip^ ). We are interested in the emission of dileptons hav- 
ing large masses and thus it is sufficient to consider emis- 
sions from the QGP phase. We also ignore the transverse 
expansion of the plasma |^] . 

As a first step, we assume that a thermally and chem- 
ically equilibrated plasma is produced in the collision at 
time tq. Assuming furthermore an isentropic expansion 
lEO. ESl one can write that 



2n^ 



1 dN 



45C(3) Ai_ dy 



(14) 



where dN/dy is the particle rapidity density for the colli- 
sion and a = 42.257r^/90 for a plasma of massless u, d and 
s quarks and gluons. We have taken the number of quark 
fiavours as w 2.5 to account for the mass of s quarks while 
evaluating a. The transverse area of the system is taken 
as A± — 7ri?^ where R± = 1.2^1^/'^ is the transverse ra- 
dius of the system for a head-on collision. We assume a 
rapid thermalization limited by tq ~ 1/3 To pl[ |. 



A±eo = / 2nre{r)dr 



(16) 



where eg is the average energy density decided by the 
temperature Tq. We note that this results in a transverse 
profile for the initial temperature as 



T{r) = To 



2 1 



r 



1/4 



(17) 



We account for the transverse profile of the jet production 
as well by introducing a factor 2(1 — r^/i?^) in expression 
(^) while performing the space-time integration 



d^x — rdr rdr drj ( 



(18) 



The limits of the r-integration are [tq, t/] for the purely 
thermal contribution. While estimating the yield of 
dileptons due to the passage of jets through the plasma 
the upper limit has to be taken as the minimum of d/c 
and Tf. Assuming a Bjorken cooling (for the chemically 
equilibrated plasma), so that T'^t is constant, we esti- 
mate Tf to be 



To 



(19) 



so that it is about 2 fm/c at SPS energies, about 3.3 fm/c 
at RHIC energies, and about 13 fm/c at LHC energies, 
when we take Tc as 160 MeV, and Tq and tq from Table |. 
Thus we note that the average value for d 



4 







C [1/GeV^] 


B [GeV] 


/3 


S 




u 


5.770x10^ 


0.2657 


6.638 


7.682 




d 


4.508x10^ 


0.2904 


6.703 


7.631 


SPS 


s 


2.760x10^ 


0.3276 


6.846 


10.696 




u 


1.147x10^ 


0.2826 


6.781 


10.949 




d 


6.028x10^ 


0.1809 


6.402 


10.942 




g 


3.180x10^ 


0.5344 


8.071 


8.964 




u 


9.113x10^ 


1.459 


7.679 


- 




d 


9.596x10^ 


1.467 


7.662 


- 


RHIC 


s 


1.038x10^ 


1.868 


8.642 


- 




u 


2.031x10^ 


1.767 


8.546 


- 




d 


2.013x10^ 


1.759 


8.566 


- 




g 


4.455x10^ 


1.7694 


8.610 


- 




u 


2.209x10* 


0.5635 


5.240 


- 




d 


2.493x10'* 


0.5522 


5.223 


- 


LHC 


s 


1.662x10^ 


0.9064 


5.548 






u 


4.581x10^ 


0.7248 


5.437 






d 


4.317x10^ 


0.7343 


5.448 






g 


1.229x10^ 


0.7717 


5.600 





TABLE II: Parameters for the minijet distribution dN/cPpr dy given in Eq. (|^) at y = for Pb+Pb at y/a^ = 17 A GeV 
(SPS), Au+Au at = 200 GeV (RHIC) and for Pb+Pb at y/s^ = 5.5 TeV (LHC). For all parameterizations K = 2.5. 

CTEQ5L plus EKS98 nuclear parton distributions are used. The range of vahdity is pr = 2 - 20 GeV for RHIC and LHC and 
2-7 GeV for SPS. 



drd (l-rVi?i) 



2 



111 3 
3 3^2 I --, -, -; 1 



3F2 



111 5 

2' 2' 2' ' 2' 



0.94 i?i 



(20) 



is smaller than the lifetime of the plasma only at LHC 
energies. In the above, the pFq's are generalized hyperge- 
ometric functions. Of course the transverse expansion of 
the plasma, which should be seen at LHC will reduce this 
lifetime ^ to some extent. Thus, on an average, the 
jet covers only a short distance in the hot plasma, before 
the medium cools down to Tc except at LHC, where 
the jets get to traverse the entire length of the system, 
while it is still hot and deconfined. We shall see that this 
leads to a large yield of dileptons from the jet-plasma 
interaction at LHC. 



C. Chemically equilibrating plasma 

Several model calculations have suggested (see e.g. 
ReL 1^) that while the plasma to be created in relativis- 
tic heavy ion collisions may attain kinetic equilibrium 
fairly quickly, the evolution to the chemical equilibrium 



proceeds only slowly through reactions of the type 
gg <-> qq and gluon multiplication and fusion gg ggg- 
Defining the quark and gluon fugacities through 

Ug = Xgflg , Uq = Xgflq , (2l) 

where Ui are the actual densities and fii are the equilib- 
rium densities for the parton species i, we can write 

TT 

^« = ^C(3)^/T3 = feiT^ . (22) 

The energy density and the pressure of the plasma are 
then given by, 

e^?,P^[a2\g + h2{\ + \g)]T\ (23) 

where 02 = Stt^/IS, &2 = lir'^Nf/m and Nf w 2.5 is the 
number of dynamical quark flavours. We further assume 
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that A„ 



Xq . Assuming a boost invariant longitudinal 



expansion, the relevant master equations can be solved 
to give 



Xg dr 



1 dXq 
Xq dr 



3 dT 
Tlh 



R- 



ai f X 



hi V A, 



where the rate constants 



i?2 w 0.24iV/a^AgTln(1.65/a^Ag 
i?3 = l.2AalT{2Xg^XlY'\ 



(24) 

(25) 
(26) 



"gl 

are taken from Ref. and include the effects of colour 
Debye screening and Landau-Pomeranchuk-Migdal sup- 
pression of the induced gluonic radiations. The initial 
conditions necessary for a numerical solution of the above 
equation are taken from the self screened parton cascade 
model 1^ and given in Table |, for the sake of complete- 



D. Jets and Drell-Yan 

We estimate the jet production and Drell-Yan produc- 
tion of dileptons in lowest order pQCD. The jet cross 
section, normalized to one nucleon pair, for the produc- 
tion of partons (quarks, antiquarks or gluons) in primary 
hard interactions between partons a and b from different 
nuclei A and B is given by [p8| 



d^PT dy 



E 

a+6— >jGt 



dXaf^{Xa)f^{x,) 

1 daa+b^lct 



■ Prey TT 



dt 



(27) 



Here and are parton distributions for the colliding 
nuclei, so that 

f^{x,Q^)^R^{x,Q^) 



7 N 



(28) 



where i?^ is the nuclear modification of the structure 
function, Z is the number of protons, N is the number 
of neutrons and = p^. The sum runs over all parton 
processes where a parton, forming a jet, can be produced 
from a and 6 with differential cross section da/dt. ^snn 
is the center of mass energy of the colliding nuclei per 
nucleon pair, so that s — Xaa;&SNN- Furthermore we have 



Xb 



XaPre 



PT^y 



PTe~ 



(29) 



(30) 
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FIG. 1: Upper part: transverse momentum distribution of 
gluons, quarks and antiquarks at j/ = for central Pb-fPb 
collisions at y^ijvjv = 17.4 GeV at SPS. Points are calculated, 
the lines are drawn as a guide to the eye. Lower part: the 
rapidity distribution for gluons, quarks and antiquarks at SPS 
for Pt = 2, 4 and 6 GeV (from top to bottom). In this, and in 
the following two figures, we have plotted the average of the 
three lightest flavours, q — {u + d+s)/3, and g — {u + d+s)/3, 
to avoid overcrowding. 



"We have used the CTEQ5L parton distributions |29f| 
and EKS98 nuclear modifications |30j. We further use 
a if-factor of 2.5 in the jet production to account for 
higher order corrections (see later). The results of our 
calculations can be found in Figs. |l|, || and ^. We also 
give a fairly accurate parameterization for the results at 
RHIC and LHC energies; 



d^PT dy 



= T, 



AA 



y=0 



da^"^ 
d^PT dy 



K 



C 



{1+Pt/B)/^ 



(31) 



for the number of quarks, antiquarks and gluons at y = 0. 
For SPS it is necessary to include a multiplicative term 
(1 - pT/8.7GeV)^ to the right hand side of Eq. (|l|) in 
order to describe the p-p-spectrum close to the kinematic 
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FIG. 2: Same as Fig. |l| for central Au+Au at ^/sjvjv = 200 
GeV at RHIC. 
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FIG. 3: Same as Fig. m for central Pb-|-Pb at -^sjvjv = 5.5 
TeV at LHC. 



limit of 8.7 GeV. Here Taa = 9A'^/8ttR]_ is the nuclear 
thickness for zero impact parameter. Values for the pa- 
rameters B, C, /3 and S can be found in Table ^ We 
note that the rapidity distributions of the partons are 
quite flat near j/ = at RHIC and LHC energies. We 
have numerically verified that assuming a flat rapidity 
spectrum at SPS energies introduces an error only of the 
order of 10-20% as the expressions favour contributions 
from rapidities near y = 0. In any case the yield at SPS 
energies is very low. 

The cross-section for the Drell-Yan process is given by, 



da 



dM^ dy 9M4 



g 

+ {q^q)] (32) 



Here the sum runs again over the quark flavours, and 



X2 = Me ^ l^/SNN- 



(33) 



These are scaled by Taa similar to Eq. ( ^l| ) to obtain 
dN/dA'P dy for the Drell-Yan production of dileptons. 



III. RESULTS AND DISCUSSION 

In Figs. ^, ^and Ijwe plot the results for thermal dilep- 
tons, dileptons from the Drell-Yan process, and the dilep- 
tons from the passage of quark jets through the plasma 
for SPS, RHIC and LHC respectively. Note that gluon 
jets will contribute only at higher order. 

At SPS energies, we recover (Fig. ^ the well known 
result that the large mass dileptons have their origin pre- 
dominantly in the Drell-Yan process. Increasing the for- 
mation time from 0.20 fm/c to 0.50 fm/c — and thus 
lowering the initial temperature by 100 MeV — drasti- 
cally alters the thermal production (from the dash-dotted 
curve to the long-dashed one) while the yield from the 
proposed jet-plasma interaction, even though essentially 
negligible, is reduced by a factor of « 2 (from the solid 
line to the long-dashed one). 

The jet-plasma interaction starts playing an interest- 
ing role at RHIC energies (Fig. ||), as now the corre- 
sponding yield is about only one third of the Drell-Yan 
contribution, and is much larger than the thermal con- 
tribution. Again lowering the initial temperature (now 
by about 150 MeV) by increasing the formation time to 
0.50 fm/c further enhances the importance of the yield 
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FIG. 4: Dilepton spectrum for Pb-f Pb at ^/s5w = 17.4 GeV 
at SPS. The long dashed curves give the resuHs when the 
formation time tq is raised to 0.50 fm/c, thus lowering the 
temperature. 



due to jet-plasma interaction. We add that this produc- 
tion is of the same order as that attributed to secondary- 
secondary quark-antiquark annihilation in a dilepton pro- 
duction calculation done using an earlier version of the 
parton cascade model [|| . 

The much larger initial temperatures likely to be at- 
tained at the LHC and the much larger (mini)jet pro- 
duction lead to an excess of large mass dileptons from 
jet-plasma interactions which can be an order of magni- 
tude larger (at M = 10 GeV) than that due to the Drell- 
Yan process. Again, reducing the initial temperature by 
raising the formation time to 0.50 fm/c reduces the jet- 
plasma yield by about a factor of 2 "while the thermal 
yield is reduced far more. We recall that at LHC energies 
several calculations (see e.g. ^ ^) have reported a 
thermal yield larger than the Drell-Yan production. We 
find that the jet-plasma interaction enhances the large 
mass dilepton production considerably. 

The equilibrating plasma scenario discussed suggests 
that, "while the quark and gluon fugacities may initially 
be less than unity, the initial temperature could be larger 
(see Table |). This results (see Fig. 0) in a situation 
"where the jet-plasma yield remains large. It is larger 
than the thermal yield, but still smaller than the Drell- 
Yan contribution at RHIC energies. 

The results at LHC energies (see Fig. ^ are partic- 
ularly interesting because the large mass yield is com- 
pletely dominated by the jet-plasma interactions, and 
this remains true even "when only jets having px > 6 
GeV are considered. 

We see that the quark jets passing through the QGP 
give rise to a large production of large mass dileptons, 
akin to the production of high energy photons through 
the same mechanism suggested earlier |l^ . They should 
be absent in pp collisions and if no QGP is formed. We 



FIG. 5: Same as Fig. |J for central Au-|-Au at -/ijvjv = 200 
GeV at RHIC. 




FIG. 6; Same as Fig. || for central Pb-|-Pb at -^/Sjviv = 5.5 
TeV at LHC. 



have noted that the results are sensitive to the initial 
conditions and the time the jet spends in the medium. 

Before summarizing, we "would like to discuss some in- 
teresting extensions of this -work. We have concentrated 
on quark jets. Gluon jets "would give rise to dileptons 
through the process qg qj*. Ho"wever, phase space 
considerations "would tend to disfavour this process "with 
respect to the channel considered in this "work. But this 
contribution may still be large, especially as the gluonic 
jets are more numerous. This could also compensate for 
the additional ag needed. However, a previous dynamical 
calculation did find this contribution to be sub-leading 
Q. The discussions in this work can be easily extended 
to diphotons Connected to that is the interesting 

question of the transverse momentum spectrum of the 
dileptons, both from the Drell-Yan process and the jet- 
plasma interaction. No K-factors were used for the elec- 
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tromagnetic processes in this first baseline investigation. 
A detailed quantitative calculation will of course require 
those, and we shall report on the results of those in- 
vestigations in a future publication |3^. We also add 
that one may in principle have a energy and px depen- 
dent K-fa,ctoT while estimating the jet-distributions. 
However, it is known that at higher energies, where the 
discussed process becomes important, the p^-dependence 
of the _fC-factor is weak , when comparing NLO con- 
tributions to the LO contributions. The precise value for 
the K-factor is thus only a constant multiplicative factor 
and results for a different value can be easily obtained. 
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FIG. 7: Dilepton yield for central Au+Au at ^/snn ~ 200 
GeV at RHIC assuming a chemically equilibrating plasma. 
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FIG. 8: Same as Fig. R for central Pb-|-Pb at ^snn ~ 5.5 
TeV at LHC. 



A more serious modification of our results is expected 



from the fact that jets passing through the QGP will 
loose energy. We note though that quark jets, the fo- 
cus of the present work, lose less energy than the gluon 
jets. Furthermore, the number of quarks in the medium 
having a given transverse momentum px drops rapidly 
as time passes, mainly due to the cooling (in the plasma) 
and also due to the free streaming of the jets. This lim- 
its the production of the dileptons mostly to the early 
times. Several authors (see e.g., Ref. Q) have discussed 
the effect of energy loss on the Drell- Yan production of 
dileptons. This can be easily incorporated in the treat- 
ment given here. 

We have not compared our results to the correlated de- 
cay of charm quarks, which provide a large contribution 
to dileptons at RHIC and LHC energies |j] , if the energy 
loss likely to be suffered by the charm quarks due to col- 
lisions and radiations is ignored. If we assume that the 
charm quarks may lose energy due to these processes , 
then these contributions will be suppressed . This de- 
bate has however not yet reached closure. A direct 
measurement could go a long way in settling this issue 
||. 

We also recall that the source of dileptons discussed in 
the present work corresponds to the so-called secondary- 
secondary interactions in the parton cascade model, dis- 
cussed in Ref. However, it was pointed out that 
the procedure adopted in Ref. Q did not account for the 
fact that the final state partons emerging from a parton- 
parton collision continue to participate in the develop- 
ment of the cascade by fragmenting and scattering while 
a virtual photon leaves the system. It will be worthwhile 
to calculate the dilepton production within a (correctly 
implemented) parton cascade model with these con- 
siderations. 

In brief, we have discussed a unique source of large 
mass dileptons in relativistic collisions of heavy nuclei. It 
arises from the passage of quark jets through the quark 
gluon plasma. The contribution is found to be largest at 
LHC energies, moderate at RHIC energies and negligible 
at SPS energies. Its detection will provide a further proof 
of the existence of the QGP as well as of multiple inter- 
actions suffered by the high energy partons eventually 
leading to jet-quenching. 
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